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Unimolecular evaporation in rotating, non-spherical atomic clusters is investigated using Phase 
Space Theory in its orbiting transition state version. The distributions of the total kinetic energy 
, release £tr and the rotational angular momentum Jr are calculated for oblate top and prolate top 

' main products with an arbitrary degree of deformation. The orientation of the angular momentum of 

the product cluster with respect to the cluster symmetry axis has also been obtained. This statistical 
approach is tested in the case of the small 8-atom Lennard- Jones cluster, for which comparison with 
7^ ■ extensive molecular dynamics simulations is presented. The role of the cluster shape has been 

systematically studied for larger, model clusters in the harmonic approximation for the vibrational 
densities of states. We find that the type of deformation (prolate vs. oblate) plays little role on 
the distributions and averages of etr and Jr except at low initial angular momentum. However, 
D ' alignment effects between the product angular momentum and the symmetry axis are found to be 

I significant, and maximum at some degree of oblateness. The effects of deformation on the rotational 

cooling and heating effects are also illustrated. 
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C/2 . I. INTRODUCTION 

^ ■ 

O • Statistical theories of unimolecular dissociation^^ provide invaluable information about binding energies, lifetimes, 
and more generally the thermodynamical state of atomic and molecular clustersj^ Besides their common use in helping 
to interpret experiments,'*'^ their quantitative relevance has been demonstrated in numerical simulationsSii with great 
accuracy. Among them, the phase space theory (PST) of Nikitin,- Light and coworkers,^ Klots,— and Chesnavich 
and Bowersii stands out as the most successfulSil in describing evaporation of weakly bound clusters for which the 
loose transition state hypothesis is satisfied. Within the PST formalism, statistical dissociation is characterized by the 
I . vibrational densities of states (VDOS) of the parent and product clusters, as well as the rotational density (RDOS) 
C ' of the products. This latter quantity depends on the interaction between the fragments, and more importantly on 
] their shapeSfii at least through their number of rotational degrees of freedom. -'^^ 

Large clusters of simple materials, such as rare gases or alkali metals, are usually found to be rather spherical. Small 
aggregates can nevertheless exhibit significant deformations at the scale of a few tens of atoms or molecules. The shape 
^ ■ of covalent clusters could be inferred from the ion mobility measurements performed by the Jarrold group In 
. ^ ■ silicon clusters, a prolate-oblate transition was evidenced around 20 atomsii^ Small carbon clusters show a competition 
, between linear chains, rings and planar shapes »^ Larger carbon clusters can be arbitrarily prolate under the nanotube 
^H ^' geometry. Deformed metal clusters have been investigated theoretically using various approacheSfi^ especially in 
O concern with the fission problem^ 

• • ' Accounting for deviations from the spherical shape may be crucial for calculating precisely the number of rovibra- 
tional states of small molecules, and the role of nonsphericity on the statistical observables relevant to unimolecular 
dissociation has been discussed previously. Chesnavich and BowersiiiS quantified the error introduced in the ro- 
^ ' tational DOS when approximating symmetric rigid rotors by spherical tops, for small molecules or radicals with 
5^ 1 moderate deformations. Bcrblinger and Schlier«i2i used a Monte Carlo integration method to compute the density of 
quantum states of the oblate molecules HJ^ or HD^ using appropriate densities or rotational states. 

Peslherbe and Hase considered evaporation in small aluminium clusters from the point of view of molecular dynamics 
(MD) simulations and PSTi In this work, deformation effects were not fully accounted for, because the centrifugal 
barrier was neglected except in the spherical top approximation for the fragments. However it is known that centrifugal 
forces, and more generally rotation may affect cluster properties quite significantly. ^Si^SiSi 
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In our previous work, an extensive study of evaporation in rotating atomic and molecular clusters was carried 
Qj^^^22j23j24 payijig a, special attention to the final angular momentum and kinetic energy released. By carefully 
calculating the vibrational and rotational DOS's pertaining to various cases of spherical clusters emitting atoms^SiSi 
linear or spherical molecules^i we could extend the conclusion of Weerasinghe and Amar^ that PST is quantitative 
in reproducing the exact numerical results. However some clear deviations could be seen in the statistical treatment 
of evaporation when applied to the smaller cluster LJg, where the main product LJ7 was approximated as a spherical 
top. This poorly satisfactory result was attributed to the non-spherical character of LJ7, known to be oblate in its 
most stable energy structure. 

Our first goal in this paper is to incorporate the effects of deformations in the PST treatment for rotating clusters. As 
a first application, atomic dissociation in LJg will be considered as an example. More generally, we will focus on large 
clusters simplified as molten, homogeneous droplets, in order to investigate the role of nonsphcricity on the observables 
relevant to cluster dissociation. Following previous effort jSiLSiiS^ it is hoped that cluster evaporation could be used as 
a probe of the shape, and possibly the phase changes, in the main product. The dynamical simulation of fragmentation 
in asymetric molecules such as those aforementioned involve can be very challenging from a computational point of 
view. Therefore it is important to incorporate and quantify the effects of deformation in the rate theories used to 
interpret gas-phase experiments on clusters. 

In the next Section, we summarize the theoretical results of phase space theory in the approximation where the 
main product after atom dissociation can be modelled as a prolate or oblate top. Concise expressions for the rotational 
densities of states, KER and final angular momentum are proposed. We also compute the distribution of the relative 
orientation between the symmetry axis and the angular momentum in the product. This property, along with its 
average value, provide an estimation of how the rotational energy is distributed in the cluster. We illustrate the specific 
case of LJg and larger LJ clusters in Sec. lIIII In particular, the roles of size, total energy and angular momentum on 
the final properties are discussed. Finally, we summarize and give some concluding remarks in Sec. II VI 



II. METHODS 



A. Phase Space Theory 

The rotational density of states (RDOS) is a key function within the framework of phase space theoryiSi^iiSiii This 
quantity depends on the kinetic energy released after dissociation, etr, sum of the translational and rotational parts. It 
also depends on the total angular momentum J. It is obtained by integration of the rotational sum of states r(e*, J^), 
which quantifies the number of rotational states available at a given angular momentum of the product cluster, 
and a kinetic energy lower than £*, its limiting value at the centrifugal barrier^ 

In PST the two fragments are treated as rigid bodies and the expressions of F only depend on their corresponding 
symmetries. Our present interest concerns the evaporation from nonspherical atomic clusters, and we have chosen to 
deal with symmetric top deformations, prolate or oblate, but with an arbitrary degree of nonsphcricity. The rotational 
constants are denoted as A and B for a prolate top cluster {X = A~ B > 0), and as B and C for an oblate top cluster 
{X — B — C > 0). The rotational energy Sr is expressed by e,. = BJ^ ± XK^, where Jr and K are the internal angular 
momenta, with -I- (resp. — ) for the prolate top (resp. oblate top) case. 

One motivation of the present work is to investigate the possible correlations between the cylindrical symmetry 
axis of the main product cluster and the angular momentum vector after fragmentation. Therefore the most general 
statistical quantity we need to consider is the probability for a dissociation event to occur with prescribed values of 
K, Jr, and etr^ 

P{K,Jr,etr;E,J)dKdJrdetrOcuj{E~EQ~etr)dKdJrdetrJ dL, (1) 

where lu is the vibrational density of states of the product cluster, Eq the energy difference between the parent and 
product, and C the available values of orbital momentum L at fixed Jr, £tr, and K. The RDOS results from summing 
the above equation over all {K, Jr) values compatible with energy and angular momentum conservation. For this we 
first need to integrate over K for rotational energies no larger than e*.^^ The expressions of these sums of states are 
given by 

2Jr =Tsie*r,Jr) if £* > ^ 

BJ^I „ , , ^ , .„ , ,T . ^ . „ ,T (2) 



A 



= Fp(£;,J^) if AJ^ >s*r> BJ^ 
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for the prolate top+atom case, and by 



2Jr — rs(£*, Jr) 



-,1/2 



A 



if e; > B J2 

rs(e;, Jr) - ro(£:, J.) if bj^ >e;> cj^ 



(3) 



for the oblate top+atom case. 

To obtain the rotational densities of states at given etr and J, the sums of rotational states must be integrated in 
the whole range C of available values of Jr and L, accounting for conservation of total energy and angular momentum. 
The general schematic integration plot is represented in Fig. ^a). Here X and Y are the two rotational constants 
of the product, such that X > Y, equal to A and B (resp. B and C) in the prolate top (resp. oblate top) case. 
The barrier height corresponding to orbital momentum L is denoted as {L) in what follows. The procedure used to 
calculate the RDOS at finite angular momentum is basically the same as used previously by Chesnavich and Bowersfii 
and by ourselves i22i24 

The intersections JriX) and JriY) have the respective values (etr/Xy^^ and (etr/i^)^^^- We introduce the limits 
L^f and LJ^ (resp. Ly and Ly) as the intersections of etr — XJr+e^{L) (resp. £tr = YJ^ +e^(L)) with L — \J — Jr\ 
(for and Ly) and with L = J + Jr (for LJ^ and Ly). At low etr, corresponding to etr < eHJ)^ and Ly are 
also given from the intersections with L = \ J — Jr\. 

Moreover we must remind that when J is not zero, the KER £tr has a minimum value e""" where integration starts. 



s™'" is such that the curve e* 



YJ^ 



This value actually depends on the smallest rotational constants, namely Y . e" 
tangent to the line L = J — Jr. 

The rotational sums of states, Eqn. Q and can be gathered into r(£*, J^) — rs(£*, J^) ± rt(e*, Jr), where 
Ts = 2Jr is the sphere+atom part, and Ft the non-spherical, symmetric top part. In the definition of F, a plus sign is 
used for the prolate top, and a minus sign for the oblate top. According to this simple partition, the RDOS can then 
be expressed as a sum of two terms: 

F(etr,J) = JJ T{e;,Jr)dJrdL 

= F,(etr, J)±rt(etr, J). (4) 

The contribution of the spherical top part to the rotational density of states has been calculated previously^SS and 
reads 



F,(etr,J) = r,-+e(£tr-S,j2)F, 



s 5 



(5) 



where 8 is the Heaviside function, and where Bs is the rotational constant equal to A (prolate top) or C (oblate top). 
The two components F" and T'^ of Fg are given by 



rs"(£tr,J) 



B. 



(J-L) 



and 



F+(£tr,J) = 2J(L. )2, 



(6) 



(7) 



respectively. 

Similarly, the non-spherical RDOS Ft(£tr, J) can be expressed as a sum of two terms, namely Ft = F^ + 0(£tr ~ 
YJ'^)T^ . These two terms cannot be fully explicited in general, due to the unknown (or complex) dependence of 
or e* with L. However, integral forms are readily available from 



Fr(£tr,J)= dL 2 

Jmax(|.I-L|,^e;/X) 



BJ^ 



1/2 



dJr 



and 



rt+(etn^) 



dL 



J+L 



L-e{6ti-X,P) Jmax(|J-L|,^6,*/X) 



BJ^ 



1/2 



dJr- 



(8) 



(9) 
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The integrals over are solved to yield 



rr(£tr,J) 



1 



dLeZ 



max \J -L\ 



' B B 



(10) 



1 f^Y 



r+(etr, J) = ^== / die; 



with the function 5* given by 



(11) 



for prolate top fragments, and by 



for oblate top fragments. 



5'(.t) = x\/ 1 — + sin ^ : 



*(x) = x\J x"^ - 1 - ln(a; + \Jx'^ - I) 



(12) 



(13) 



B. Kinetic energy release distribution 

From the density of rotational states, the kinetic energy release distribution can be easily calculated from the usual 
relationship 



P(etr;-B, J) cx / T{et,,J)uj{E~EQ~et,)detr 
in which ld corresponds to the vibrational density of states of the product cluster. 



(14) 



C. Product angular momentum distribution 

We keep the notations introduced previously for the rotational constants Bs, X and F, such that {Bg^X^Y) = 
{A,A,B) for the prolate top fragment, and {Bs,X,Y) = {C,B,C) for the oblate top fragment. We also note 
Lx — (e^)^^(£tr — XJ^), where (e^)^^ is the reciprocal function of e^. 

The probability distribution of Jr, P{Jr', E, J) is obtained from integration of Eq. Q over K and etr- This quantity 
has been investigated in details in our past workJi A similar approach can be pursued here, after using the fact that 
the RDOS is the sum of the two contributions Fg and ±Ft. Hence the distribution P can be written as P = Pg ± ^t, 
where the plus (resp. minus) sign is used for the prolate (resp. oblate) top fragment. The expressions for the spherical 
and symmetric top contributions are given by 



and 



PtiJr;E, J) 



Ps{Jr;E,J) = 2Jr e[LBA£tr,Jr)-\J-Jr\] 

xuj{E- Eo - etr)[niin(J + Jr,isJ - \J - Jr\]deti-, 



(15) 



2 

7x 



E-Eo 



e[LYiStr;Jr)-\J-Jr\] 
^inin( J+ Ji-,Li') 



xu;{E - Eo~ etr)rf£tr 



let-BJ^l^^^dL. 



<(|J-J,|,Lx) 



(16) 
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D. cos 6' distribution 

The relative orientation of the angular momentum of the product cluster with respect to its symmetry axis can be 
calculated by integrating Eq. (QJ over all variables except K. At given Etn Jr and L, the number of available K states 
depends on the values of e* = etr — e^-^) and J^. In the prolate top case, there are two states available if e* > AJ"^. 
If AJ"^ > £r — BJ^, there are also two states available provided that \K\ < i^max = [(e* — BJ^)/X\^/^, and zero 
otherwise. The relative orientation cos9 — \K\/Jr must then satisfy 

e;>{Acos^9 + Bsm'^0)J^. (17) 

The same ideas lead to the following condition for the oblate top case: 

< > {Ccos'^e + Bsin^e)J^. (18) 

The integration over the {L, Jr) plane is now restricted, as indicated in Fig.QJb), to include the above conditions. At 
a given cos 6, the intersection Jr{9) with the L = axis is given by Jr{0) — [etr/ (B sin^ 9 + Bg cos^ 6)Y^^ in general. 
This quantity continuously sweeps the [Jr(X), Jr{Y)] interval. 

Finally we get the probability density that the relative orientation between the angular momentum vector and the 
revolution axes has the value cos 9 within d cos 9 as 

P{cos9;E,J) (X / Lj{E - Eo - etr)detr X \ 2dLdJr+ 

26 [e; - {B sin^ 9 + B, cos^ 9)J^] dLdJr \ , (19) 

X .J^>e^>Y J'^ J 

valid for both types of deformations. These distributions, as well as the preceding probabilities of £tr or J^, must be 
calculated numerically in general. Ingredients other than the rotational densities, most importantly the vibrational 
densities of states and the centrifugal energies s^{L), are obtained from the same techniques described in our previous 
work,^^ 

However, the general shape of the distributions can be inferred already at this stage. In the prolate top situation, 
Bg > B and P{cos9) is maximal when cos0 — 0, and minimal when cos9 =^ 1. The reverse holds for the oblate 
top case. Therefore, P{cos9) is a decreasing (resp. increasing) function of cos6' for prolate (resp. oblate) fragments. 
This means that the angular momentum gets more likely aligned with the principal axis with largest inertia. Thus it 
simply reflects the greater rotational stability of rigid bodies when their angular velocity is loweriSS 

III. RESULTS AND DISCUSSION 

The clusters we are interested in are Lennard- Jones clusters characterized by an arbitrary degree of deformation. 
The simplest model is that of an ellipsoid, whose axes ratio can take any real value. At a given (large) number of 
constituants n, the reasonable assumption of constant density implies that the volume of the cluster is constant, which 
allows us to define the equivalent spherical radius Rg as Rg = Ron^^^, where Rq plays the role of a lattice parameter. 
In practice, Rq is taken to reproduce the rotational constant of the icosahedral 13-atom cluster. 

The ellipsoidal shape of the cluster is continuously varied by keeping its volume constant. By denoting a and b the 
long parallel and short perpendicular axes of the ellipsoid, respectively, we quantify deformation in the cluster using 
'J — (b ~ a)/{b + a). Equating i?f to ab^ yields 

^=l + (a/i?.)3/2- 

As the cluster becomes increasingly prolate, 7 decreases and tends to —1, which corresponds to a linear system. For 
oblate systems, 7 increases toward +1, which corresponds to a planar system. The rotational constants corresponding 
to the deformation 7 are given by 

and 



B = 2Bo 



1 + 7 



4/3 + 



(22) 
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A. LJg — > LJt+LJ 

The dissociation of LJg has been considered first since it provides a good candidate to test the proposed formahsm, 
the product cluster LJ7 being not perfectly spherical in its ground state geometry. For this system, the anharmonic 
VDOS and the effective dissociation potential were obtained from parallel tempering Monte Carlo and Wang-Landau 
simulations, respectively. The technical details can be found in our previous work^ 

To quantify the real extent of deformation in the vibrationally excited LJ7 cluster, we have analysed the thermal 
evolution of its rotational constants. In Fig |21 the average constants are plotted versus T. At T=0, the minimum 
energy configuration is the pentagonal bipyramid, an exactly oblate system for which A = B ^ 0.17, C « 0.11. These 
rotational constants yield a deformation index 7 « 0.3, or a/Rs = 0.7 with Rs = 0.55. At higher temperatures, 
the cluster deviates more and more from the oblate shape, eventually becoming prolate. In order to get meaningful 
statistics of the dissociation process using standard molecular dynamics, we had to thermalize the parent cluster at 
rather high temperatures, typically above the melting point (T « 0.2). The rotational constants used in the PST 
analysis were taken at this precise temperature, yielding A « 0.18 and B = C ~ 0.11 LJ units. These values give an 
equivalent deformation index 7 — —0.3 in the ellipsoidal picture, or a/ Rs — 1.5. Interestingly, the double icosahedron 
LJig is more spherical than LJ7 (7 « —0.18), and remains so even above its melting temperature. 

In Table^we have reported the values of (etr) and (Jr) obtained by MD simulations and by the PST descriptions 
in various rigid body approximations for LJ7. These results correspond to E/n — 1.2 and J = 0, 1, 2 and 3 (one LJ 
unit of angular momentum approximately equals 337i for argon). The data for J = is obtained using simplifications 
to Eqn. H14l - I19|l . which arise due to the new restraint on the orbital momentum, namely L = —J^iii^'i^'^^^ 

Under these aforementioned conditions, increasing cluster deformation (7 < or 7 > 0) mainly increases the 
product angular momentum, and marginally reduces the KER. The prolate top approximation performs significantly 
better than the simple spherical assumption in reproducing the MD results, but this is also true for the oblate top 
approximation. 

The complete distribution of Jr contains more information than the average value alone. This distribution is 
plotted in Fig. O for the initial angular momentum J = 3. A good agreement between MD and PST is obtained 
when the deformation of the product cluster is taken into account. By comparison, the distribution in the spherical 
approximation underestimates Jr- More interestingly, the PST/oblate calculation looks very close to the prolate result, 
suggesting that the type of deformation does not alter significantly the rotational distribution after evaporation. The 
KER distribution, not plotted here, shows a similar behavior. However, it should be kept in mind here that we are 
dealing with only moderate deviations from sphericity. 

More insight is gained from the distribution of cos 6, represented in Fig. 0| for the same conditions of total energy 
and angular momentum. The distribution is obviously uniform in the spherical case, as there is no privileged axis. 
As expected from our previous analysis, the prolate description of the product cluster favors low cos9. This is in 
agreement with MD, but differs from the PST result in the oblate approximation. 

Increasing angular momentum in the parent cluster brings the final rotation axis closer to the major principal axis, 
as shown in Fig. [S| This emphasizes the need for properly describing deformations in clusters having a significant or 
even moderate rotational motion. The values of J considered in this figure are small enough to produce rotational 
heating However high angular momenta, which result in rotational cooling, also display a similar behavior. 
Comparing the values of (cos 9) at J = and J = 3 shows that the deviation from random orientation (for which 
(cos 61) — 1/2) is half due to the extent of deformation, and half due to angular momentum itself. From this figure, 
but also from Fig. ^ alignment between the rotation axis and the main symmetry axis is about twice more efficient 
for prolate tops than for oblate tops. Such an effect is expected to be also magnified by the extent of deformation. 

B. Size effects in model clusters 

We now discuss the general effects of deformation on the statistical properties of larger clusters after evaporation. 
Our interest here is not to reproduce MD or experimental results, for this reason we have used an approximate model, 
where the n-atom LJ cluster is treated as a continuous medium with constant density, and the vibrational density of 
states is assumed to be harmonic. Such a model allows us to vary arbitrarily the degree of deformation 7 in the whole 
range —1 < 7 < 1. The interaction potential is also simplified as — 4n/r^, which yields a barrier energy proportional 
to L^. As deformation of the cluster increases, the physical extent of the cluster also increases and we expect the 
1/r^ approximation to fail eventually. A more appropriate potential such as l/(r — tq)^ would not be relevant for the 
present, qualitative discussion, even though it would provide a more robust ground for eventual comparisons. 

Three parent cluster sizes n+1 have been selected, which correspond to n ~ 50, 100, and 200, respectively. To cover 
a vast range of situations, two total energies {E/n = 0.9 and 1.2 LJ units per atom) and two angular momenta ( J = 
and J = 20 LJ units) have been considered, thus providing four different physical conditions for fragmentation. The 
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effects of deformation on the kinetic energy released and on the final angular momentum are shown in Figs. and 
respectively, for the four present situations. From a general point of view, smaller clusters show more pronounced 
deviations of the average KER, product angular momentum. This is easily understood by noticing that the rotational 
constants, hence all the parameters which characterize the deformation such as A, decrease with size very fast (oc -nT^^^) 
with respect to, e.g. the centrifugal barrier energy (cx n~^/^). In other terms, the specific symmetric top contributions 
7t to the RDOS become less important with increasing size. 

It is useful to analyse first the J = case, since all the integrals involved are one-dimensional. For convenience, we 
assume in the following discussion that the main product is a prolate top. In this case the RDOS is exactly given for 
J — {) and a Cq/t^ interaction potential 



r(£tr,J = o) = 



J, 



2 J y* dij 



+ / 2 



etr - a - BJ^ 
A 



1/2 



dJr, 



(23) 



where aJ^ — e^{Jr) is the centrifugal energy at the barrier. Here a does not depend on 7 but only on n and Cq. 
and are solution of £tr = AJ^ + aJ^ and £tr = BJ^ + aJ^, respectively. When a = 0, this integral can be exactly 
solved as 



r(£tr, J = 0) = 



£tr 



(24) 



r is a linear function of £tr independently of the rotational constants, therefore the average KER does not change 
with deformation, because of the normalization of the integrals involved in the average. This shows that any variation 
in the average energy released is due to the centrifugal barrier only. The Klots models-'^^ approximate the rotational 
densities r(£tr,low J) as a power law in £tr, the exponent increasing proportionally with the number of rotational 
degrees of freedom. Since this number does not depend on the possible symmetry of the products, there are no 
deformation effects on the evaporation statistics within the Klots model. 

We have not found simple expressions for the RDOS by Taylor expanding the general formula of Eq. H23fl in a, 
because the result still contained complex functions of the rotational constants. Instead we consider the Jr observable. 
The general expression for the probability distribution of at J = 0, for prolate top product, readsS^ 

P{Jr;J = 0) ^ 2Jr uj{E - Eo - etr)detr 

JAJ^+aJ^ 



BJ^+aJ^ 



£tr 



A 



Uj{E - Eq- £tr)d£tr 



(25) 



After some algebra, the lowest order in expansion in |7|"'^^^ of this expression is found to be 

pE—Eo 

P{Jr;J = 0) ~ 2Jr / uj {E - Eo - etr)detr 

JBoJ^+aJ^ 

+ 2V2B^^^\-/\'^/^J*uj'iE -Eo- BoJ^ - aJ^), 



(26) 

However, because the 



where uj' — dio/dE. As I7I increases, the weight of the term increases as well as {Jr 
correction has the power 3/2, the deviation with respect to the sphere remains small. 

A similar treatment can also be carried ou in the case of oblate top products. However, the result differs because 
the smaller rotational constants changes from B ~ Bo{l + 47/3) for 7 < to Bs — Bo{l — 27/3) for 7 > 0: 



P(J^;J = 0) ~ 2Jr 



E — Eq 



uj{E - Eo- sti-)detv 



^-fjME -Eo- BoJl - aJl), 



(27) 



The corrective term is still positive, and favors higher values of a-t increasing 7. However, this term grows linearly 
with 7, therefore the correction is more important for oblate tops than for prolate tops. 
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The increase in product angular momentum, as well as the non equivalence of the two deformations close to the 
spherical shape, are indeed seen in Figs.EIa,b), and these effects are also observed at nonzero initial angular momenta, 
Figs. [3c, d). Obviously, other parameters to quantify the extent of deformation might affect the shape of the curves 
{Jr)il) or (etr)(7)- For instance, if 7 was defined from the inertia momenta or rotational constants directly, (J^) 
would vary differently around zero, because the inertia momenta vary quadratically with the axes lengths a and b. 

Increasing I7I tends to favor larger product angular momenta, which is a consequence of smaller rotational constants, 
the same rotational energy being achieved through higher Jr . In turn, the centrifugal barrier increases and the kinetic 
energy released decreases. This explains the correlation between (etr) and (Jr) seen in Figs. [HI and [T] Larger total 
energies decrease the relative importance of the rotational contribution, as well as the effects of nonsphericity. This 
is precisely what we observe in Figs.EJb) and^lb). Conversely, if now we increase J keeping the total energy fixed, 
more kinetic energy is allowed in the rotational modes, which usually yields amplified deviations with respect to the 
spherical reference for both the KER and the product angular momentum. 

At finite angular momentum J, it is much harder to provide simple explanations for the behavior of the statistical 
quantities, because new parameters come to play an important role. Increasing J generally has a much more dramatic 
effect than increasing E. The minimum value for the KER, e™'", has a component which grows linearly with the 
smallest rotational constant of the product. ^^■■^'^ At large deformations, either prolate or oblate, e}."'" always takes 
values smaller than for the spherical reference. Therefore low values of the KER are further statistically favored. This 
effect may be larger than the small increase due to the centrifugal energies seen in Fig. |HIa). Second, because high 
deformations are characterized with at least one very small rotational constant, the upper limits Jr(Y) of the product 
angular momentum in the [L, Jr) integration range can reach higher values. Thus the integration range extends to 
larger J^, and the average value of (J^) increases. 

The competition between these various effects, as well as the different functional forms for the RDOS, provide a 
rather rich variety of behaviors for (etr) and (Jr), as shown in Figs.Eland0 The exact numerical calculations displayed 
in Figs. EIc,d) and[7Ic,d) confirm that much larger deviations are seen with respect to the spherical reference when 
J is nonzero. In particular, the marked decreases in the KER and the correlated increases in the product angular 
momentum are of larger magnitude than for J = 0. A notable result is that the most important differences between 
the two types of deformation (oblate and prolate) are exclusively seen at very low total momentum and for large sizes. 
These differences are related to variations in small rotational constants. Above some specific total energy and/or some 
total angular momentum, the effects of prolate or oblate deformations become rather similar in shape and magnitude, 
except very close to 7 = 0. 

Another interesting feature of these nonspherical systems is their ability to adopt a preferential direction for the 
rotational excitation following evaporation, as measured with respect to the cylindrical axis. The average cosine of 
the angle between the product angular momentum and the symmetry axis, as obtained from Eq. (|I9II . is shown in 
Fig. IHlfor the same mechanical conditions as in Figs. El and [7| versus the extent of deformation 7. Again, the effects 
monotonically decrease with increasing cluster size. 

From the previous discussion we already know that (cos 6*) is positive (resp. negative) for oblate (resp. prolate) 
deformations. Hence, the variations of this quantity seen in Fig. (Hjare not surprising, at least for small deformations. 
A simple, first-order perturbative expansion in the equation H19|l can be used as a check of the behavior of (cos 9) 
close to 7 ~ 0. However, for arbitrary total energies and angular momenta, the extent of alignment (cos0) always 
shows the same qualitative variations with respect to 7. For oblate top products, dissociation induces rotation around 
a long axis preferentially, and the effect is more significant at smaller total energies or larger initial angular momenta. 
Prolate top products, on the other hand, display a non-monotonic behavior, concomitant with the variations of the 
large rotational constant, Eq. 1)221) . which has a maximum when 7 = {V^ — l)/(-\/2 + 1)- The fact that alignment is 
hindered at very large oblate deformations is a simple consequence that the two rotational constants B and C tend to 
0, hence the second integral over L and Jr in Eq. (|I9|I goes to 0, and the system becomes equivalent to a sphere again 
from the point of view of its orientation. Therefore, it is practically impossible to reach a perfect statistical alignment 
for oblately deformed systems. Such a situation is never met with prolate deformed clusters, since one rotational 
constant always keeps increasing for large deformations. The shape of the curves in Figs. |Htc,d) also reveals that 
alignment can be quite effective, even at moderate deformations. For a 50-atom cluster, J ~ 20 is close to the typical 
thermal angular momentum near the melting point, and a deformation of about 20% is enough to induce statistical 
alignment beyond 25%. 

Finally, following our previous work^^i we focus on the influence of nonsphericity on the rotational cooling and 
heating effects. The same simple model for LJ clusters was used in the harmonic approximation for the VDOS and 
the C/r^ radial potential for the dissociation energy. We consider here more realistic conditions for dissociation, 
namely thermal distributions for the energy and angular momentum of the parent cluster, as well as a possible extra 
energy shift to model an initial, brief excitation (from photoabsorption or collision). To provide ground for comparison, 
these distributions are assumed to be Boltzmann-like and identical for all geometries. The final distributions were 
calculated by taking a full account of the prolate or oblate character of the main product cluster. 



9 



In spherical clusters, we have seen that a purely thermal distribution generally leads to rotational cooling, i.e. the 
distribution of angular momentum of the product cluster is shifted to lower values^ On the other hand, an initially 
cold cluster submitted to a sudden vibrational excitation preferentially exhibits rotational heatingiSiS& In Fig. 1^1 the 
final angular momentum distributions are compared for the three types of products. Even though the product clusters 
are significantly deformed, the general rotational behavior remains similar to that of the spherical top. As seen from 
Fig.EIa), purely thermal evaporation induces rotational cooling, and this effect is attenuated for nonspherical clusters. 
Deformation also influences the final distribution itself, which differs significantly from the Boltzmann behavior for 
prolate systems. These results are consistent with our previous analysis, as thermal evaporation corresponds more 
likely to the cases (c) and (d) of Figs. EHHl with substantial initial J. At I7I = 0.5, Fig. |7fd) shows that (J^) increases 
by about 5% for oblate top products, slightly less for prolate top products, with respect to the spherical case. These 
changes (including their magnitude) are reflected on the distributions in Fig. (SJa). 

Rotational heating is produced in a different way, by adding a sharp extra vibrational excitation to an initially cold 
cluster. Deformed clusters also exhibit rotational heating, as shown in Fig. I^Jb). For both prolate and oblate top 
products, deformation further heats the final rotational motion, but the two deformations behave dissimilarly, the 
effect being much stronger for the oblate top product. This is understood by considering again Fig.Efb), which shares 
the mechanical conditions of high total energies but low angular momenta. The increase in (J,.) is quite important 
for oblate tops, although it is barely noticeable for prolate tops at I7I — 0.5. 

Deformation most often induces higher rotational excitations upon evaporation, this explains why it reduces rota- 
tional cooling, but enhances rotational heating. 

IV. CONCLUSION 

The present work was aimed at providing a quantitative study of the effects of cluster shape on its evaporation 
statistics. Within the framework of phase space theory, the distributions of kinetic energy released and final angular 
momenta were obtained assuming the fragments could be described as an atom plus a prolate or an oblate top as the 
main product. 

These theoretical tools were challenged on the example of the unimolecular dissociation of the 8-atom Lcnnard- 
Jones cluster, chosen for its significantly oblate shape. The small disagreement reported earlier^^ between PST and 
molecular dynamics simulations for the distribution of final angular momentum was solved after accounting for the 
nonspherical character of LJ7. 

Using simple models for the vibrational density of states and the dissociation potential, we systematically investi- 
gated the changes induced by deformation. We found very significant deviations with respect to the spherical case 
for both the final angular momentum and kinetic energy, especially for nonzero initial momenta. In such cases, the 
relative variations of the energetic properties are larger than those of angular momenta by up to even one order of 
magnitude. 

We could interpret these changes as originating partly from the specific centrifugal barrier in our present imple- 
mentation of PST, but also from the new energy constraints arising with lower rotational constants. These effects are 
not included in simpler statistical rate theories of unimolecular dissociation such as the Klots models, not mentioning 
the Engelking-Weisskopf^-'^° or RRK approaches?^ 

We also obtained information about the relative orientation of the product angular momentum, by considering 
restricted integration of the rotational densities of states. Alignment toward the long axes is generally observed for 
all deformations, in agreement with mechanical arguments. This alignment was seen to be very effficient already at 
moderate initial angular momenta, and for moderate deformations. 

Finally, we discussed the influence of deformation on the rotational cooling and heating effects resulting from specific 
excitations. While cooling is attenuated for nonspherical clusters, heating is amplifled. These behaviors are attributed 
to the decrease of one rotational constant, which favors larger angular momenta in the product cluster. 

The examples considered here were obtained on model clusters, but we expect the present results to be relevant 
for more realistic systems. Charged rare-gas clusters, which display very prolate shapes below 15 atoms due to the 
presence of a linear ionic trimer, have been the subject of intense experimental and theoretical activity from the point 
of view of their dissociation patterns. ^"^ 

Polycyclic aromatic hydrocarbon (PAH) molecules are strongly oblate and their rotation was considered within the 
astrophysical context by Rouan and co-workersi^ Dehydrogenation of these molecules, in particular, has received a 
special attention among several experimental and theoretical groups^^*^^*^ This problem could be investigated with 
the same methods used here. 

In these two examples, a realistic atomistic description requires models much more involved than simple pairwise 
potentials. This practical limitation prevents one fully relying on brute-force simulations, thus making statistical 
approaches most useful in the future. 
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0.89 
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0.89 
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TABLE I: Average kinetic energy released (etr) and product angular momentum (Jr) after evaporation of the LJs cluster at total 
excitation energy E/n = 1.2 and several total angular momenta. The values obtained from molecular dynamics simulations 
arc compared to the PST predictions in the spherical (7 = 0), prolate top (7 = —0.3) and oblate top (7 = 0.3) approximations 
for the product cluster LJ7. 
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FIG. 1: Range of integration in the {L,Jr) plane for atom+symetric top fragments, (a) Complete integration range, X >Y 
being the two rotational constants, (b) Restricted integration range when the orientation of Jr with respect to the symmetry 
axis is restricted to have a fixed angle 9. 




0.1 0.2 0.3 0.4 

Temperature 



FIG. 2: Tliornially averaged rotational constants of the LJ7 cluster versus temperature, for initially nonrotating (solid lines) 
and rotating (J = 3, dashed lines) systems. 
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FIG. 3: Distribution of final angular momentum after dissociation of LJg, from diflferent approximations (sphere-atom, prolate- 
atom and oblate-atom) for the rotational DOS in the PST calculation. 




I I I I I I I I I I I I I I I I I I I L 

0.2 0.4 0.6 0.8 
COS 



FIG. 4: Distribution of cos 6 in the dissociation of L,J8 at J=3 and E/n=\.2. The results of MD simulations are compared to 
PST calculations for different approximations of the LJ7 main product. 
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FIG. 5: Average angle (cos^) between the symmetry axis and the final angular momentum versus initial J in the dissociation 
of LJs with E/n=l.2. 
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FIG. 6: Average KER (etr) versus deformation 7 in the evaporation of model LJ„+i clusters under several conditions, (a) 
E/n = 1.2 and J = 0; (b) E/n = 1.2 and J = 20; (c) E/n = 0.9 and J = 0; and (d) E/n = 0.9 and J = 20. The results are 
shown for n = 50 (solid lines), n = 100 (dashed lines), and n = 200 (dotted lines). 




FIG. 7: Average product angular momentum (Jr) versus deformation 7 in the evaporation of model LJ„+i clusters under the 
same conditions as in Fig.ini The same graphs conventions are also used. 
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FIG. 9: Initial and final distributions of angular momentum after evaporation from a 51-atom model LJ cluster, the main 
product being assumed spherical (7 = 0.0), prolate (7 = —0.5), or oblate (7 = 0.5). (a) Hot thermal distributions of energy 
and angular momentum (T = 0.5); (b) cold thermal distributions (T = 0.05) and extra vibrational excitation {AE = 30). 



